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Starting from the rigorous excitation equation, the propagation of waves through a 2D waveguide 
with the periodically corrugated finite-length insert is examined in detail. The corrugation profile is 
chosen to obey the property that its amplitude is small as compared to the waveguide width, whereas 
the sharpness of the asperities is arbitrarily large. With the aid of the method of mode separation, 
which was developed earlier for inhomogeneous-in-bulk waveguide systems [Waves Random Media 
10, 395 (2000)], the corrugated segment of the waveguide is shown to serve as the effective scattering 
barrier whose width is coincident with the length of the insert and the average height is controlled by 
the sharpness of boundary asperities. Due to this barrier, the mode spectrum of the waveguide can 
be substantially rarefied and adjusted so as to reduce the number of extended modes to the value 
arbitrarily less than that in the absence of corrugation (up to zero), without changing considerably 
the waveguide average width. 

PACS numbers: 41.20.-q, 84.40.Az, 78.67.Lt, 89.90.+n 



I. INTRODUCTION 

The scattering through waveguides with surface rough- 
ness is a ubiquitous phenomenon which occurs on differ- 
ent length and time scales both in the natural environ- 
ment and in different artificial systems |Q, ||. Whereas 
for environmental waveguides the problems related to 
wave scattering at surfaces of random nature are of major 
importance, for man-made guiding systems the particu- 
lar attention is called by the scattering produced by the 
boundaries with periodical corrugation. 

Periodic structures, in general, have received much at- 
tention in different areas of physics through the variety 
of applications of their peculiar spectra. In particular, 
waveguides with periodic inserts have for a long time 
been employed as slow-wave structures in vacuum and 
accelerator physics photonic (Bragg) crystals are 

widely used as frequency selective narrowband resonators 
and filters in terahertz and optical ranges ; multilay- 
ered films with periodic variation of the refractive index 
serve as effective mirrors in the VUV and X ranges of EM 
radiation ||, |^. The importance of periodic structures 
for numerous applications and the diversity of physical 
effects originating in them have stimulated the develop- 
ment of a number of theoretical methods for their anal- 
ysis. The detailed review of the methods developed up 
to the mid fifties of the 20-th century can be found in 
Ref. 0. 

For periodic structures, apart from the foremost ques- 
tion of how the forbidden gaps in their spectra originate. 



there arises a number of other intriguing questions re- 
garding their physical properties, especially for those sys- 
tems whose inhomogeneities are much less in size than the 
operating wavelength. Until recently it was believed that 
inhomogeneities of this sort cannot crucially affect the 
wave propagation. For this reason the studies on struc- 
tures with sub- wavelength periodic variations of their pa- 
rameters have for long remained beyond the main line of 
research until the previous decade, when anomalous wave 
transmission through sub-wavelength hole arrays in op- 
tically thick films was detected Q [ill- In our recent 
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paper |13 it was also discovered that the small-scale dis- 
tortion of wave resonator boundaries dramatically affects 
their spectra. Specifically, it was found that modal con- 
tent of rough-side resonators was mainly controlled by 
the sharpness of boundary asperities rather than by their 
height. Even for the asperities of sub-wavelength nature 
the spectrum may be adjusted in such a way that at any 
desirable frequency range the quasioptical resonator can 
become even a single-frequency one. 

The results of Ref. were obtained for closed sys- 
tems only, being related to their spectral properties. It 
would be, however, of great interest to find out if small 
boundary asperities could produce similar radical effects 
on the properties of open systems, in particular, on the 
wave transmission in waveguide structures. Yet, in try- 
ing to answer this question one faces the problem which 
has not so far been adequately resolved. Specifically, any 
self-consistent theory of a waveguide system with arbi- 
trarily corrugated side walls is relied unavoidably on the 
results related to plane-wave scattering by such bound- 
aries. For boundaries whose corrugation is random such 
theories do exist for a long time. Yet they are all known 
to allow obtaining the tractable and practically useful re- 
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suits only in the case of smooth roughness, i.e., for the 
asperities whose mean-square tangent of slope is small as 
compared with unity For cavity resonators, this 

restriction was overcome in Ref. [|l3| , but for systems of 
general type, including waveguides and different surfaces 
subject to radiation, this is yet to be done. 

In studying the wave scattering by distorted bound- 
aries the most challenging task is, in general, related to 
boundaries' random character. The scattering by the 
surfaces whose corrugation is deterministic, in particu- 
lar, periodic in tangent coordinate seems to be, at first 
glance, a less intriguing problem. Normally, the peri- 
odicity as such is considered to be the main feature of 
such systems, whereas the the specific details like the 
amplitude and the sharpness of corrugation receive less 
attention. By now, there has been a number of works on 
wave propagation in periodically corrugated waveguides 
(see, e.g., Refs. [p5|-p0[ and references therein). Many 
peculiarities of their spectra are studied in more detail. 
But in all of the aforesaid papers the particular effective- 
medium models were used. In spite of the fact that these 
models allow to take adequately into account the peri- 
odicity property, they were not capable of allowing for 
other substantial features of the wave scattering by un- 
even surfaces. 

Meanwhile, in Refs. 
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through randomly rough single-mode quantum waveg- 
uides it was revealed that the particular role in the 
transport and spectral properties of such systems was 
played by the so-called "slope" , or "gradient" , scattering 
mechanism, in which as the dominant factor for surface- 
roughncss-induced scattering of propagating waves serves 
the degree of asperity smoothness. The same scattering 
mechanism was found to influence radically the spectra 
of cavity resonators with sharply rough walls . As far 
as the principal difference between random and deter- 
ministic corrugation lies in the difference in correlation 
lengths, the results of Refs. |2^, ^ suggest quite 
strong grounds to believe that the gradients of boundary 
roughness will have a great impact on the properties of 
periodically corrugated waveguides as well. 

To verify this conjecture, in this study we solve the 
problem of wave transport through the waveguide with 
periodically corrugated finite-length segment whose cor- 
rugation is thought of as being small in amplitude but 
arbitrarily large in sharpness. Through the application 
of special operator technique of transverse-mode separa- 
tion the transport of each of the transverse mode in the 
waveguide with corrugated insert is shown to be equiva- 
lent to quantum-particle transport through the individ- 
ual sharply modulated potential barrier whose average 
magnitude is governed by the degree of asperity sharp- 
ness. The small-scale nature of the asperities makes it 
possible to operate in the lowest Bloch band only, which 
allows to concentrate the attention on the sharpness- 
induced scattering leaving the band-forming Bragg scat- 
tering beyond consideration. Within the suggested ap- 
proach it is revealed that the slight-in-amplitudc modu- 



lation of the waveguide boundaries, taken to be not ex- 
tremely smooth, affects significantly the waveguide mode 
structure. On sharpening the asperities the number of 
propagating modes reduces substantially, the fast modes 
turning sequentially, first, into the slow and then into 
evanescent modes. In such a way the multimode waveg- 
uide can be made single-mode or even evanescent-mode 
one without noticeable variation of its local width. 



II. THE PROBLEM FORMULATION 

We will consider the infinitely long two-dimensional 
waveguide with straight parallel boundaries, whose side 
walls at some finite segment of length L are subjected 
to periodical corrugation (see Fig. |l|). The width of 
the entire waveguide will be specified by the function 
w{z) — wa + 2£^{z)9{L/2 — \z\), where wq is the width of 
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FIG. 1: (Color online) Two-dimensional 
waveguide with periodically corrugated insert 
(segment z e [-L/2,L/2]). 

the straight segments / and ///, ^(z) is the periodic func- 
tion with period d, amplitude a and zero mean value. To 
simplify further numerical calculations, but without sub- 
stantial loss of generality, we will describe the corrugation 
by function ^(z) = cr cos(27rz/c?), taking the amplitude 
small as compared to waveguide average width, 

<T<$:wo . (1) 

Consider the experimental setup when in a real 3D 
waveguide the TE type waves are excited, namely, the 
TEon modes with indices corresponding to field varia- 
tions along X and y axes, respectively. In this case, all the 
electromagnetic field components can be found with the 
aid of one (x) component of vector potential A(y, z) that 
obeys the homogeneous wave equation. The waveguide 
side boundaries will be considered as ideally conduct- 
ing, which for s-polarized waves corresponds to Dirich- 
let boundary conditions at surfaces y — ±'w{z)/2 while 
on the straight a;-coordinate boundaries (not shown in 
Fig. if they are taken into consideration, the New- 
mann conditions should be met. As far as the conditions 
are concerned at the open waveguide ends, we will discuss 
them below, after the main equation for the sought-for 
field will be written down in the mode representation. 

To solve the excitation problem for the waveguide we 
will seek for the Green function which obeys the following 
equation 

{A + k^)G{r,r')^S{r~r') , (2) 
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where A is the two-dimensional Laplace operator. As thus expressing the Green function as a double series 
a first step, we Fourier transform equation (^) with re- 
spect to the locally (in the coordinate z) complete set of 
eigenfunctions of the transverse Laplacian, viz. 



G(r,r')= ^ {y,z;n|G„„,(^,z')l2/',^';"') • (4) 
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For mode coefficients G„„'(z, z') the following set of cou- 
pled equations results from Eq. 
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where ^ — (ttji/wo)^ is the longitudinal propaga- 
tion number for the n-th waveguide mode in the non- 
corrugated outward parts of the waveguide, Vn{z) and 
Unm{z) are the effective "potentials" (in general, of op- 
erator nature) which arise in the wave equation due to 
the boundary corrugation. 
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These potentials are concentrated across the finite inter- 
val, so we will regard them as the sources for intra- and 
intermode scattering of outward waveguide modes, re- 
spectively. The numerical coefficients in potential ( |6b| ) 
have the form 
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Equation (g|) must be supplemented with boundary 
conditions (BC)'s along axis z. Since the waveguide is 
an open-ended system, as the BC's at z — > ±oo we will 
use the radiation conditions. To formulate them ana- 
lytically for a non-homogeneous waveguide system is a 
challenging task in the general case. However, in view of 
the perturbation potentials in Eq. (^ being identically 
equal to zero in the regions \z\ > L/2, the Green function 
mode structure in these regions reduces to its diagonal 
components. So, the legible boundary conditions may 
be a priori formulated for these components only. The 
next section will show that this is sufficient for determin- 
ing the entire Green function since all the non-diagonal 



components of its mode matrix can be lineally expressed 
in terms of diagonal elements, for which we obtain the 
whole set of uncoupled equations. 



III. MODE SEPARATION IN THE WAVEGUIDE 
WITH NONUNIFORM SEGMENT 

The set of equations (|5|) may be solved for mode com- 
ponents of the Green function using the special operator 
technique applicable for potentials Vn{z) and Unm{z) of 
quite an arbitrary form. The technique was first elabo- 
rated in Ref. |24j for the open systems of waveguide con- 
figuration and then expanded to open-ended and closed 
resonator-type systems disordered in the bulk as well 
as on the surface ||l3, 26 . With regard to the waveguide- 
type systems, the method reduces (schematically) to the 
following action sequence. At the first stage, by putting 
n 7^ n' in the set of equations (||) one should resolve this 
set with respect to the off-diagonal components of mode 
matrix ||Gmn||, thereby expressing them, by means of lin- 
ear operation, through the diagonal elements whose mode 
indices are coincident with the right-hand (column) index 
of the required off-diagonal element. 



Gmn{z,z') = J dziKmn{z,Zi)Gnn{zi,z') [m ^ n) . 

(8) 

In Eq. (^, the integration runs over the axis z region 
where potentials (^ are not equal to zero. The kernel 
of the integral operator in (^) can be found from the 
equation of Lippmann-Schwinger type, 



y^ran{z, z') = Rmn{z , z')V^ j dzi Rrak{z, Zl ) Kfe„ (zi , z') 

(9) 



where we introduce the notation 

R„„(z,z') = G(^)(z,z')C/™„(z') 



(10) 



for matrix elements of some operator R whose domain of 
definition is two-dimensional space M = {z, to} which in- 
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eludes the coordinate axis (z) and the whole set of mode 
indices (m). This operator apparently accounts for scat- 
tering between all the different transverse modes, so in 
what follows we will refer to it as the mode-mixing oper- 
ator. 

Function Gm"^ (z, z') in expression ([l0|), which we iden- 
tify as the trial mode Green function, is the solution of 
equation 



Q2 

— + kl~Vn{z) 



G'y\z,z')^5{z-z') , (11) 



which differs from Eq. (|^) by the lack of intermodc poten- 
tials (to unify the notations in the subsequent formulas 
we re-designate the mode variable of the Green function 
in Eq. ( pi] ) in comparison with Eq. (|lO|), viz. m — > n). 

Upon solving equation ( |ll|) and considering the oper- 
ator R matrix elements as already known functions, one 
can substitute off-diagonal elements of the Green matrix 
into equation (^, thereby obtaining closed equations for 
the diagonal elements of this matrix (for Vn's), 
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Gnn{z,z')=S{z^z') . (12) 



In equation (|T^), if compared with Eq. (|TT|) for the trial 
Green function, the additional operator potential Tn has 
appeared, whose expression has the form 



(13) 



Here, t( is the intermode operator potential which is spec- 
ified in M by matrix elements 



{z,n\U\z',n') = C/„„/(z)(5(z - z') 



(14) 



P„ is the projection operator, whose action reduces to the 
assignment of given value n to the nearest mode index of 
any operator standing next to it, no matter if it is to the 
left or to the right. Two embracing projectors in Eq. ([l^ ) 
have appeared due to the restriction of summation in 
Eqs. (|^) and ^ by mode indices k ^ n. The operators 
standing between the projectors act in the subspace M„ £ 
M which includes the coordinate axis z and the set of all 
mode indices other than the separate index n. The role of 
the projectors in operator potential Eq. ( [l^ ) is to reduce 
the action of this operator on the function standing to the 
right of it solely to integration over z, without summation 
over mode indices. 

From the functional structure of potential ( p^ ) one can 
infer that in equation ( |l2| ) for the diagonal propagator of 
a given transverse mode this potential accounts actually 
for the intermode scattering. In contrast to potential 
Vn{z), the potential Tn, in the general case, acts as the 
non-local (in z) operator whose characteristic scale is de- 
termined by the spatial extent of the trial Green func- 
tions and the intermode potentials, both entering the T- 
potcntial through the mode-mixing operator. 

As regards the boundary conditions to equation (p^, 
the n-th waveguide mode, whose dynamics is governed 



by this equation, is proven to be effectively separated 
from other modes, which are present in Eq. (^^ only im- 
plicitly, as the intermediate scattering modes "hidden" 
inside the T-potential. If Eq. (^2|) is recast in the form of 

the integral Dyson equation, where G^\z, z') stands as 
the unperturbed sought-for solution, it will be clear that 
open boundary conditions for the exact intramode Green 
function coincide with those satisfied by the correspond- 
ing trial function. The conditions for the latter function 
will be presented in the next section. 

Thus, taking account of the relationship (^) between 
diagonal and off-diagonal elements of the mode Green 
matrix, which can be recast in the operator form 



Gn 



Prn{t-R)-'RP„Gr. 



(15) 



with Gmn and G„„ being the operators in the coordinate 
variable z, one can assert that the solution of the entire 
set of uncoupled equations (p2|) (for Vri's) completely de- 
termines the sought-for Green function of the waveguide 
with the inhomogeneous insert. 

Upon deriving equation (^2]) and explicit expressions 
for the potentials entering it, we in fact have performed 
the effective "mode separation" in the original two- 
dimensional Helmholtz equation. This operation sub- 
stantially facilitates the solution of the problem of cal- 
culating the fields in the originally two-dimensional in- 
homogeneous waveguide since there exists a number of 
advanced mathematical methods to solve wave equations 
in the dimension one. 

However, it should be admitted that at the expense 
of the achieved simplification, in the master equation, 
whose role from now on will be played by equation (12), 



we have obtained the additional (to the original Vn{z)Jm- 
tramodc potential with quite non-trivial functional struc- 
ture. The effective potential of the analogous type is 
familiar in the quantum scattering theory ||2^, 28 , be- 



ing referred to as T-matrix. However, the conventionally 
defined T-matrix is known to be an extremely singular 
mathematical object, whose calculation is normally car- 
ried out in the lowest order of perturbation theory. In our 
calculation technique the T-potcntial in Eq. (12) has no 
singularities. In Ref. |2^ it was shown that if the effective 
potentials in Eq. (|^) are random functions, which break 
the spatial symmetry of the system under consideration, 
the regularity of potential (13) is ensured by splitting the 



entire set of effective potentials into the subsets of intra- 
and intermode ones. But the system we consider in this 
study is subject to regular rather than random pertur- 
bations, so it is not a priori clear if the arguments of 
Ref. are appropriate for this case. 

At the same time, the system we deal with is the open 
system and thus non-Hermitian, which normally is suf- 
ficient for the inverse operator in Eq. (p^ ) to be nonsin- 
gular. The proof of this fact can be found in fundamen- 
tal papers by Feshbach on nuclear reaction theory (see 
also Rcfs. where this issue is discussed in detail 

with regard to one-dimensional systems subject to final- 
support external potentials). Yet, to additionally protect 
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oneself against the degeneracy property of the full wave 
operator in Eq. (||) one can always supply the "energy" 
k'^ with additional infinitesimal positive-imaginary term 
which accounts for the dissipation and certainly ensures 
the lack of deg eneracy, thus rendering the inverse opera- 
tor in Eq. ( [l3[) properly defined. 



IV. CALCULATION OF THE INTRAMODE 
PROPAGATORS 



As far as the sought-for Green function (^) is com- 
pletely determined by the diagonal components of its 
mode matrix, from here on we will regard the solution 
of Eq. as the main object of our study. Under 

fixed mode index n this equation can be thought of as 
being uncoupled provided that all trial Green functions 
Gm \z, z') with mode indices m ^ n are beforehand de- 
termined. The next section is devoted just to the solu- 
tion of this auxiliary problem, namely, to finding the trial 
functions. 



A. The trial Green function 

Function G^\z,z') solves one-dimensional equation 
( pi] ) on the condition that both ends of the waveguide 
arc open. This condition, referred to in electrodynam- 
ics as Sommerfeld's radiation condition |3^, can be 
formulated in the form of two equalities connecting the 
sought-for function and its first derivative, viz. 



dz 



= 



(16) 



Being so posed, the problem of finding the trial Green 
function pertains to the class of boundary-value problems 
of open type, whose solution is difficult to obtain even for 
potentials of relatively simple form. 

Fortunately, in the case of one spatial dimension the 
boundary-value Green function problem reduces to the 
pair of problems of causal type (see, e. g., Ref. |Q). This 
can be done through representation 

Gl^)(z,z') =[W„(z')]"' [^+(z|n)V'_(z»e(z-z') 

+ ^+(z»^_(z|n)e(z'-z)] , (17) 

where functions ^± (z|n) are the solutions of two auxiliary 
Cauchy problems for homogeneous equation (pH), whose 
boundary conditions are specified on solely one ("plus" 
or "minus" , respectively) end of the interval where Green 
function is being sought for; yV„(z) is the Wronskian of 
those solutions, which, in concordance with Liouville's 
formula for equation (^ij) with potential ( |6a| ) does 
not depend on the coordinate z. 

In this paper we choose the potential Vn{z) to be a pe- 
riodic function in the range z G [— L/2, L/2]. Basically, 
this enables one to explore the solution of Eq. (|l^) tak- 
ing advantage of Floquet theory p6i. Yet, the structure 



of the right-hand side (r.h.s.) of Eq. ( pa] ) is excessively 
complicated for this purpose, so we will make some sim- 
plifications allowing for the reduction of Eq. (11) to the 
standard Mathieu equation. 

We will first regard geometric parameters of corruga- 
tion as obeying the inequalities 



cr, d < Wo 



(18) 



At the same time, the sharpness of corrugation will be 
considered as arbitrary. By separating in the potential 
( |6a| ) its mean value and the term that oscillates about 
zero, Vn{z) — Vn{z) + l^Vn{z)j under conditions (JTsj) we 
have 



Vn{z) = 

AK(^) 



2^2 



^2 



27r2 



cos 



A-KZ 



(19a) 
(19b) 



where parameter S = a/d is introduced, which specifies 
the degree of the asperity sharpness. The sum of poten- 
tials (^9|) is shown schematically in Fig. ||, making us to 
expect that the corrugation of the finite waveguide por- 
tion would result in the effects substantially analogous to 
those arising in the quantum particle transport through 
periodically modulated potential barrier. The particular 




FIG. 2: The effective potential barrier in equation (|ll| ) 
arisen due to the waveguide corrugation over the finite 
interval L. 

feature of our problem is that the average height of the ef- 
fective potential barrier, being proportional to the square 
of parameter S, appears to be governed by the sharpness 
of the corrugation profile. Below we give a brief account 
of the calculation procedure for causal functions '0±(z|n) 
and, correspondingly, for Green function (p7[). 

Consider, for definiteness, function i/'+i^T"-); function 
'!/'-(z|n) can be found analogously. By going over to 
the dimensionless coordinate t = 27rz/d the equation for 
the sought-for function reduces to the canonical Mathieu 
equation 



V'+(i|n) -f [an - 2g„cos(2t)]V+(t|n) = 



(20) 
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Here, prime signs stand for the derivatives over <-variable 
and the notations are used 
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(21b) 



Bearing in mind that at the "plus" end of the waveg- 
uide fimction ip^(z\n) obeys radiation condition similar 
to Eq. (|l^) we will seek this function in the form 

^+iz\n) = < An(Pn{z) + Bn^Pni-z) , \z\ < L/2 , 

[t„e*'="(^-^/2) ^ z>L/2 

(22a) 

where 
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FIG. 3: Stability areas (grey) for equation Eq. (gO|) so- 
lutions. The shading approximately covers the region 
corresponding to inequalities (25). 



"5(^1") 



(22b) 



is the standard Floquet-Bloch solution of equation ( |20| ) 
on the entire (infinite) axis z 5 , Tji and tji are 

the constants subject to determination from the match- 
ing conditions at the junction points between corrugated 
and straight waveguide portions. Bloc h w avenumber 
>f„ — 2TT^n/d in the exponential in Eq. (22b) is related 
to characteristic number ^„ of equation (20), being the 
function of parameters a„ and g„ ||37| . Functions (/3„ (±2;) 



in the second line of Eq. (22a) arc the linearly indepen- 
dent solutions of the wave equation on the interval of 
waveguide corrugation. Function g{z\n) from Eq. (221;) 
is, in the general case, representable as a series. 



g{z\n) 



Jo 



(n) 



E 



(23) 



whose coefficients are specified by the parameters of 
equation ( pO| ) and obey the following set of recurrence 
relations, 



[an 
(m 



(/A 
0, 



> + 2™)^]/^'^ 
±1, ±2, ...) 



/•(") 



(24) 



In the present study, we focus our attention on the 
effect produced on the spectral and dynamical proper- 
ties of the waveguide by the gradient scattering mecha- 
nism per se, without making an emphasis on spectrum 
multi-band structure associated with the periodic nature 
of waveguide corrugation. Therefore, in what follows we 
will restrict ourselves by the values of parameters a„ and 
g„ that fall into the shaded region in Fig. ^ corresponding 
to the couple of inequalities. 



a„ < 1 



l9n 



(25) 



The region includes the extended states of only first sta- 
bility band of solutions to equation (]20|), as well as the 
evanescent states with negative mode energies. In terms 



of physical parameters of the system at hand, the in- 
equalities (|2^) suggest the relative smallness of the corru- 
gation amplitude and the smallness of its period against 
the wave length, i. e., the fulfillment of inequalities 



Nr. 



A 



2-K 



> d 



(26) 



where Nc is the number of extended modes in the cor- 
rugated region. Under conditions (|2^), the terms of se- 
ries ( ^3| ) decrease rapidly with growing |m|, thus permit- 
ting us to restrict the summation by three terms only, 
with indices m = 0,±1. In this case the characteristic 
number of equation ( pO| ) differs insignificantly from its 
value at (jn = 0, so in the estimations one can always put 

/^n ~ \fan- 

The number of extended trial modes in the corru- 
gated segment of the waveguide is determined by the 
condition o f posi tiveness of the entire mode "energy", 
k\ — y,i(z), and therefore is largely dependent on 
the corrugation profile sharpness. 



A^c(S) 



/V2 . 



2~? 



l + (27r2/3)S2 



(27) 



In Eq. (|^), square brackets denote an integer part of 
the enclosed number, Nco = [kwo/Tr] is the number of 
extended modes in the corrugation-free waveguide parts. 
In Fig. H, the diagrams of A^c(S) are presented for several 
specific values of the average waveguide width. It is eas- 
ily seen that on sharpening the asperity profile the num- 
ber of extended modes in the corrugated region decreases 
quite rapidly, reaching zero at some critical value of the 
sharpness parameter S, namely, at Scr = Nco/V2. The 
wavenumber >f„ in Eq. (22b) is real- valued at 5 < Scr, 
while at S > Scr the corrugated segment of the waveg- 
uide becomes practically "non-transparent" , in the sense 
that all modes within it turn into the evanescent ones, 
each decreasing exponentially over the scale of the order 
of |>f„|-^ 
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(28c) 
(28d) 



In Eqs. (28c) and ( P8d| ) we have introduced the notations 



.g\L/2\n) 



g{L/2\n) 
= g{L/2\n) . 



(29) 



FIG. 4: The number of extended modes in the corru- 
gated waveguide segment versus the sharpness of cor- 
rugation profile. The diagrams correspond to Nco = 
5, 15, and 30. 



By moving along axis z in Fig. |l| from region III in the 
negative direction and matching the fields sequentially at 
boundaries between regions III-II and II-I we obtain the 
set of equations for coefficients entering Eq. (22a). The 
solution of this set reads, 



^ _ (l - Ul) exp (iKnL) 
1 — T^^j exp (2ix-„L) 



Tin 



1 — exp (2iKnL^ 
1 — 7?.^ exp {2ixnL) 



(28a) 
(28b) 



Note that in deriving Eqs. (^8|) the waveguide side bound- 
aries were assumed to be continuous at points z = ±L/2, 
but not necessarily smooth. It can be easily verified 
that coefficients t„ and r„, Eqs. (28a) and (281), re- 
late with one another via equality + Vn]"^ = 1, which 



suggests the validity of their interpretation in Eq. (22a) 
as the mode transmission and refiection coefficients, re- 
spectively. 

The calculation scheme for function ip-{z\n) is com- 
pletely analogous to that for ip^(z\n), differing only by 
the inversion of axis z. By performing the calculations 
and by inserting the obtained expressions for 'ijj±{z\n) 
into formula (p^), the trial Green function in the corru- 
gated waveguide segment can be expressed as a sum of 
four terms each of which is proportional to the product of 
Bloch eigenfunctions running in opposite directions, viz. 



Gi^)(z,z') = t;i^)e*-"(^-^')g(z|n)5(-z'|n) + f;(^)c-^'^"(^-^')g(-z|n)g(zV) 

+ c;(^)e-"(^+^')g(z|n)g(z» + CJ^^^e— "(^+^')g(-z|n)5(-z» . (30) 



Factors C/j^''' in Eq. (^) have the form 



Mn) _ _ 
»11 ^ 



. Qn 



9{z - z') + 7^2 e2'-"^0(z' - z) 



iz' -z) + nle^^'-^eiz - z') 



Mn) _ Mn) _ . Qn ^ 

and the notations are used 



7^r 

Qn = 



1 



(31a) 

(31b) 
(31c) 

(32a) 
(32b) 



Parameter TZn in Eq. ( p2aD obviously resembles the 
well-known expression for the reflection coefficient from 
the potential step, which is familiar from standard quan- 
tum mechanics textbooks (see, e.g., Ref. [Q). It only 
differs from the latter coefficient by the slight modifica- 
tion (by substitution >Cn — > /3n) related to periodicity of 
the potential on the one side of the step in the case un- 
der examination. In our particular problem, the stepwise 
jump of the effective potential results from the gradient 
renormalization of the mode energy while crossing the 
border between regions I or III of the waveguide and the 
corrugated region II (see Fig. |^) wherein the mean value 
of the potential differs from zero. 

The quantity Q„, which enters the Green function 
terms as a common multiplier, is of interest in itself, so 
it is worthwhile to examine it in more detail. This factor 
has a resonant character and is modulo maximal when 
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the Bragg condition is met 



It: 



(33) 



which means that exactly half-integer number of the n'th 
mode wave lengths are accommodated on the length of 
the corrugated segment. In the plots of |Q„| versus fre- 
quency (see Fig. bearing in mind that through the 
Green function the mode density of states is expressed, 
one can see the well-known van Hove singularities, which 
are related to the mode cutting-off in going across the 
boundary between straight and corrugated regions of the 
waveguide. 
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FIG. 5: (Color online) The plots of the quantities \Qn\, 
|7?.„P and versus the frequency measured in units 
of ttc/wq. The points of TZn slumps in the upper plot 
correspond to the opening of the next waveguide chan- 
nel. Geometrical parameters of the corrugated segment 
are chosen as follows: S — 3.33, a — 0.04 wg, L = 0.3 



For the purposes of comparison, in Fig. [s] the depen- 
dencies on dimensionless frequency ZU = (^{wo/nc) are 
shown for three specific modes (n = 1,2,3) of (a) reso- 
nance factor Qn, which determines the intensity of the 
field of the corresponding mode "accumulated" in the 
corrugated segment, (b) the modulo squared reflection 
coefficient from the boundary between corrugated and 
straight segments of the waveguide (TZn), and (c) of the 
modulo squared coefficient t„ of the mode transmission 
through the corrugated region. It is seen that at the 
points of appearance of the next propagating mode the 
mode density of states (proportional to Q„) has obvi- 
ous singularity. To the left of the corresponding cut- 
off frequency (Un^), the coefficient of mode reflection 
from the effective potential discontinuity (TZn), in view 
of the mode assuming the evanescent character, is of the 
modulo one, whereas the transmission coefficient turns 
to zero. To the right of a;„ the TZn modulus rapidly 
drops down, whereas the transmission factor abruptly 
rises, being at the same time subject to resonance oscil- 
lations | |4l[ |. The ideal transparency of the waveguide in 
the oscillation maxima (|i„| = 1) corresponds to the ap- 
pearance inside the corrugated segment of quasi-bound 
states, which arise due to the interference of waves par- 
tially reflected from the interfaces between waveguide 
segments with a different mode structure. Thus, the cor- 
rugated segment of the waveguide plays actually the role 
of an open-ended resonator in which, due to its openness, 
the unlimited energy accumulation even at frequencies so 
close to the resonance ones is forbidden. 



B. The exact mode propagators 

The above obtained trial Green function may be 
thought of BjS cl starting point for obtaining the exact 
mode propagator G'„„(z,z') using the perturbation the- 
ory with regard to the intermode scattering. Such an 
interpretation of function G^n^\z, z') is corroborated by 
the fact that both of the Green functions are connected 
with one another by Lippmaim-Schwinger equation, see 
Ref. . The exact mode propagator may differ substan- 
tially from the trial one, since in Eq. ( p^ ) the effective 
potential Tn has appeared whose effect on the waveguide 
spectrum is not necessarily small. From the functional 
structure of this potential one can infer that within the 
framework of our theory the intermode scattering reduces 
to the trial mode intermixing, which is implemented by 
operator R in Eq. (p^. 

To analyze the intermode scattering in the arbitrary 
order of intermode potentials, with regard for the po- 
tential (^3|) complex structure, is a rather difficult task. 
This can be accomplished efficiently only in the limiting 
cases of weak and strong trial mode intermixing. The 
intensity of the intermode scattering may be quantified 
by the operator R norm, depending on whose value we 
may scrutinize both of the above indicated limits by ex- 
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panding the T-potential in an appropriate manner. 

In Appendix |^ the operator R norm is estimated at 
an arbitrary value of the sharpness parameter 5. It 
is shown that with the chosen waveguide parameters, 
specifically, if the length of the corrugated segment ex- 
ceeds considerably both the wavelength and corruga- 
tion period (L ^ k~^,d), while the parameter S does 
not exceed its critical value Scr, the operator R norm 
appears to have, for the most part, very large values, 
||R|| ^ 1. It only assumes small values, as compared to 
unity, when either the corrugation is extremely smooth, 
so that S ^ {kL)^^{d/ L), or if parameter S has ovcrcrit- 
ical values, where the corrugated segment of the waveg- 
uide becomes the evanescent-mode one. However, in the 
dominant range of parameters, including the case where 
the corrugation is rather smooth (S ^ 1), the inequal- 
ity holds ||Rj| S> 1, as is exemplified in Fig. || where the 
numerical results are presented for the operator R norm 
calculated for the particular set of waveguide parameters. 



be transformed in the following way. 




FIG. 6: The operator R norm versus the degree of the 
corrugation sharpness. The comb-shaped structure of 
the graph is related to the sequential closing of the ex- 
tended modes, in compliance with Eq. (p7|). The waveg- 
uide parameters are as follows: a — O.OBwq, L — IOwq, 
N,o = 15. 

At first glance, the fact that the trial modes become 
substantially intermixed, even if the corrugation is rather 
smooth, should imply that they cannot serve as a good 
approximation for determining the exact mode states. 
Yet, in the case of the large operator R norm one can 
carry out the following formal manipulations with equa- 
tion ( |l^ ) . By considering this equation as the coordinate 
matrix element of operator equation 



(34) 



where is the operator whose coordinate matrix el- 
ement is the Green function (|30|), one can expand the 
inverse operator in Eq. (|l3|) into a series not in the oper- 
ator R itself but rather in the its inverse, R~^. The ex- 
pression located between the projectors in Eq. (nSf) may 



R-Ul- R) 



-u{i - R-^y 



(35) 



Taking into account that operator U is of zero-diagonal 
nature in the mode variables, one can reduce the poten- 
tial Eq. (f3|) to the form 



Tn 



l + (R-i) 



(36) 



By substituting it to equation (p^, we arrive at the 
following operator expression for intramode propagator 

Grin (^! ^ ) ■ 



Gr. 



i(R-i) 

2 ^ ^ nn 



or in the coordinate representation 



(37) 



(38) 



Based upon this result it may be concluded that if 
the corrugation-induced intermodc scattering becomes 
strong, in terms of the trial mode intermixing, the precise 
waveguide spectrum returns with parametric accuracy 
to the form specified by the initial trial Green function, 
which allows for the intramode scattering only. Consid- 
ering the intermode potentials in Eq. (^) for ||R|| 3> 1 
results in the twofold reduction of the amplitude of the 
wave excited by a given source in comparison to the am- 
plitude of the wave we have chosen as the trial one. 

Interestingly, although the intermode scattering in the 
system of our study is, chiefly, strong this should not have 
a strong impact upon the spectra of signals transmitted 
across the corrugated waveguide segment. As seen from 
relation (|l5|), at large operator R norms for propagator 
Gmn{z, z') with m 7^ n the operator (in z coordinate) 
expression is valid 



Gn 



R- 



. Gn 



(39) 



From Eq. (|39|) it follows that if the waveguide is excited 
by the mode which is not of beyond-cutoff nature for 
corrugated segment, at the output of this segment the 
mode should be basically registered of the same mode 
index. The other modes should have small amplitudes at 
the exit of the corrugated region, in spite of the high-level 
intermode scattering within it. Although this conclusion 
seems to be unexpected at first sight, it is fully consistent 
with the data obtained, e. g., by Vellekoop and Mosk |49 
in their studies of wave propagation through waveguides 
with a variable level of inhomogcneity, as well as with 
qualitative explanation of their results suggested later on 
in Ref. l50|l. 
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V. CORRUGATION-INDUCED 
RAREFICATION OF WAVEGUIDE SPECTRUM 

The above consideration suggests that regardless of 
the corrugation-induced intermode scattering being weak 
or strong, the spectrum of the corrugated waveguide is 
mainly determined by the set of trial mode Green func- 
tions, i. c., it can be with good accuracy approximated by 
the entire set of Bloch wave numbers, For corrugated 
waveguide walls, the dependence of these wavenumbers 
on frequency lo bears the same qualitative character as 
for the ideally straight-wall waveguide. Yet, there is a 
number of essential differences between waveguides with 
straight and corrugated boundaries. One of the differ- 
ences we have revealed is that the cutoff frequencies of 
the corrugated waveguide exceed the corresponding fre- 
quencies of the waveguide with straight walls. This can 




5 10 15 20 

FIG. 7: (Color online) The Brillouin diagrams for cor- 
rugated (dashed lines) and the corresponding non- 
corrugated waveguide (solid lines). The corrugation 
parameters are: a = 0.05 wq, S = 0.25. 

be seen from Fig. |^, where Brillouin diagrams resulting 
from our formulas are presented for three different modes 
of the corrugated and the corresponding non-corrugated 
waveguide. 

Besides, even small-amplitude corrugation of the 
waveguide walls has a noticeable effect upon the den- 
sity of its cut-off frequency spectrum. Without cor- 
rugation, the set of these frequencies is equidistant in 
mode numbers, w^'^'' ~ T:nc/wo. If corrugation is present, 
the equidistance is violated and the interval between the 
neighboring cut-off frequencies increases with a growth in 
the corrugation sharpness. This is clearly evident from 
Fig. 1^, where the relative interval between neighboring 
cut-off frequencies is depicted, namely. 



where fin is the n-th mode cut-off frequency for the cor- 
rugated waveguide. One can see from the graphs that the 
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FIG. 8: (Color online) The relative interval between 
neighboring cut-off frequencies for the corrugated 
waveguide against the corrugation sharpness. Param- 
eter (T is fixed and equal to 0.05 wq- 



spectrum of modes propagating in the waveguide with 
corrugated walls is considerably rarefied with an increase 
in the corrugation sharpness, as compared to the simi- 
lar spectrum of the waveguide of the same average width 
whose side boundaries are ideally smooth. 

The result we have obtained is in conformity with the 
general claim regarding the normal frequencies of con- 
servative two-dimensional systems , namely, that the 
addition of couplings into a dynamical system can only 
enlarge the intervals between its normal frequencies. In 
our particular case, this statement turned out to be valid 
for infinite-dimensional dynamical system (^, where the 
role of coupling is played by the level of corrugation. Note 
that the effect of the cut-off frequency shift to a higher- 
frequency region was already noted in Ref. Yet, it 
was not appropriately examined therein and no general 
rules were established. 

One more intriguing property of the corrugated waveg- 
uide is closely related to the cut-off frequency growth, 
namely, a decrease in the number of modes that are al- 
lowed to pass through the waveguide, as the sharpness of 
corrugation increases with its amplitude being kept con- 
stant. As can be seen from Eq. (|27|), a decrease in the 
number of extended modes with increasing the sharpness 
parameter S may be quite significant, up to the waveg- 
uide being entirely "locked out" , while the number of 
modes being allowed to propagate in the non-corrugated 
regions is reasonably large. This property may be used 
for creating the waveguide filters, as it permits quite 
significant spectral modifications of signals transmitted 
through artificially corrugated waveguides. 
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VI. THE ENERGY TRANSPORT IN THE WAVEGUIDE WITH CORRUGATED SEGMENT 



The Green function obtained in Section IV, apart from the spectral properties of the waveguide with corrugated 



insert, aUows for direct calculations of its conducting properties, in particular, of power transmission. The power 
radiated by the plane antenna positioned at cross-section Sz' and then transferred through cross-section Sz at the 
point with longitudinal coordinate z, being averaged over oscillation period is given by formula 



w(z)/2 w{z)/2 



P{z,z') ^ -CklmJ ds J ds'G{r,r')^G*{r,r') ^ - {Cll)klm J dy J d?/'G(r, r')^G*(r, r') 



(41) 



-■w{z)/2 ~w{z)/2 



Here, C is the dimensional constant, whose presence relates to the fact that the dimensionality of Green function from 
equation does not meet that of the electric field from Maxwell's equations; is the waveguide dimension along 
axis X, the asterisk denotes complex conjugation. By substituting Green function in the form of double series (^) into 
Eq. ( pl| ) we arrive at the expression 



P{z,z') = -{Cll)kIm J2 Gmniz,z')—G*„,,,iz,z') 

m.n— 1 



ci 



oo oo 



- k Im / / 
w{z) ^-^ 

\ ^ n=l i,fc=l 



2ik 



i+k 



(42) 



In the context of our problem the second term in the r.h.s. of Eq. ( p^ ) may be omitted. This is substantiated, 
first, by the observation that this term does not contain the products of the diagonal mode Green functions. As far 
as the off-diagonal ones are concerned, we neglect them in view of equality (^9[). The additional argument to neglect 
the second term in the r.h.s. of Eq. (^) is the presence of factor w'{z)/w{z), which is exactly zero if the receiver is 
positioned outside the corrugated region. With allowance for the result ( p8| ) we thus arrive at the following formula 
for the power transmitted across the corrugated segment of the waveguide, 



iv). 



z, z 



d 

— ( 

dz' 



z = L/2 
z' = -~L/2 



4c 



E 



ikn 



(43) 



To make a detailed calculations of this powe r, one should substitute either expression ( |30| ) for the trial Green 
function or the transmission coefficient Eq. (28a) into formula (^). In figs. H and 10, the power transferred through 
the corrugated segment of the waveguide is depicted as a function of the normalized frequency and the sharpness 
parameter, S. The stepwise form of both of the graphs is related to the addition (or, correspondingly, cutting-off) 
the effective (gradicnt-rcnormalized) extended modes in the waveguide spectrum, in accordance with formula (|27[). 




6 8 10 12 

FIG. 9: (Color online) Frequency dependence of the 
power transmitted across the corrugated waveguide 
segment (ct = 0.02wo, L = l.bwo, S = 1). 




FIG. 10: (Color online) The transmitted power 
against the corrugation sharpness (cr = 0.02u'o, 
L = 2wo, oj = 6ttc/wo). 
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However, whereas the steps in Fig. |9| are the well-known physical fact^the analogous stepwise behavior of the trans- 
mitted power as a function of corrugation sharpness, shown in Fig. |lO|, is predicted therein for the first time. The 
oscillations seen on the plateaus in both of the graphs, ^ and |o[ are related to Bragg resonances in the transmission 
of the marginal waveguide mode. 



VII. SUMMARY AND DISCUSSION 

To conclude, we formulate the main results obtained 
in the present paper and give some comments on the 
exploited calculation technique. We have developed the 
theory of wave transmission through the waveguide con- 
taining the finite segment with periodically corrugated 
side boundaries. Such an inhomogcneous insert in the 
initially unfigured waveguide can have a considerable im- 
pact on its spectral properties. The central role in the 
modification of the entire waveguide mode structure is 
played by the sharpness of the insert corrugation. Even 
though the boundary asperities might be of very small 
amplitude, the number of extended modes in the waveg- 
uide is substantially reduced when the corrugation is 
sharpened, reaching zero at some critical value of the 
governing sharpness parameter. 

To solve the problem we have applied the method of 
oscillation mode separation, which was developed earlier 
with regard to waveguide-like systems with bulk inho- 
mogcneities. In the present study we have adjusted this 
method to waveguides with inhomogcneous boundaries 
by using thereto the local transverse-mode representa- 
tion, which makes it possible to present surface inhomo- 
geneity as the bulk one. 

The separation of waveguide modes, which is basically 
equivalent to the separation of variables in the original 
wave equation, was carried out through the specific oper- 
ator procedure essentially equivalent to the construction 
of the perturbation series with regard to the intensity 
of the intermode scattering. In developing this theory, 
we started from the solution to the auxiliary problem, in 
which the intermode scattering is formally disregarded 
while the intramode scattering is taken into account in 
full. The application of this procedure has enabled us to 
find that the intermode scattering induced by the cor- 
rugated waveguide segment is normally very strong, re- 
sulting actually in the mode system homogenization. In 
these conditions, the spectrum of sharply uneven waveg- 
uide is restored to the state were the intermode scatter- 
ing is entirely absent while the intramode scattering is 
accounted for. The actually large strength of the inter- 
mode scattering reveals itself in the twofold reduction of 
the field amplitude in comparison with the case were this 
scattering is disregarded arbitrarily. 

By solving the set of decoupled equations derived using 
our method, we succeeded in showing that the corrugated 
segment of the waveguide may be considered (by analogy 
with quantum particle transport) as the modulated po- 
tential barrier, whose width is coincident with the actual 
length of the segment whereas the height is prescribed by 
the level of corrugation sharpness. Inside this sharpness- 



controlled barrier the field is transferred in a resonant 
manner by means of the quasi-steady states arising due 
to the multiple wave reflection from the barrier ends. The 
"external" waveguide modes falling upon the barrier are 
mutually re-scattered while passing through it, but at the 
exit point their arrangement in mode indices is basically 
restored. As the barrier increases with growing the corru- 
gation sharpness, the propagating modes are successively 
cut off, though the asperity height is insufficient for this 
purpose. The mode cutting-off leads to rarefication of the 
entire (combined) waveguide spectrum, so that at some 
critical level of the asperity sharpness the waveguide be- 
comes the evancsccnt-modc one. 

Another important result of the application of our 
method is that by using it one can analytically trace the 
overall dynamics of the propagation constants against the 
sharpness of boundary corrugation. Although functional 
dependence of the longitudinal wavenumber on the fre- 
quency for both the straight and the corrugated waveg- 
uide remains the same, in the presence of corrugation 
the cut-off frequencies line up non-equidistantly. The 
shift of these frequencies to the higher-frequency range 
grows considerably both with a decrease in the corruga- 
tion period and with a mode number increase. The dif- 
ference between neighboring cut-off frequencies increases 
as well, so that the waveguide spectrum becomes signifi- 
cantly more sparse. 

The non-equidistance of the cut-off frequency set 
may be utilized in practice to achieve the single- 
frequency mode in microwave amplifiers, in developing 
the waveguide-based filters for complex signals, as well 
as for narrow-band generators with Bragg reflectors. In 
such generators, the amplification curve normally has fi- 
nite width. Therefore, in order to prevent the given mode 
interaction with higher modes it is necessary that they do 
not fall into the amplification band. At the same time, 
the operating mode frequency, at which the output-input 
ratio reaches a maximal value, is situated at a relative dis- 
tance of approximately 1/4 from the amplification curve 
lower bound. So, it is quite desirable for the cut-off fre- 
quencies to be shifted disproportionately to mode num- 
bers. 



Appendix A: Estimation of the mode-mixing 
operator norm 

According to the standard definition of the operator 
norm MS, H, we have to calculate the quantity 



(R^, R^) 

sup ^7 r-^ 



(Al) 
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where trial functions tp{z) must belong to the functional 
space where operator R is active. The domain of def- 
inition of this operator, as shown in Sec. Ill, is the 



coordinate- mode space, so the result of its action on some 
scalar function (in particular, on the single component of 
a vector function) is a vector. 

The support of potentials (||) is the segment C : \z\ < 
L/2. Therefore, it is exactly this particular segment that 
we will consider a s th e support for the entire set of trial 
functions in Eq. ( |AlD . The finite-support functions are 
Fourier-expandable, so in the capacity of trial functions 
in the definition Eq. ( |AlD it is reasonable to use the set of 
exponential functions ipi^) = exp(igz). The compactness 
of the support for the latter functions may be ensured 
by choosing the integration contour for Fourier integrals 
to be not the exacly real axis, q £ M., but rather the 
horizontal line in the plane of the complex variable (q) 
properly offset to the upper or lower half-plane, c. g., q Cz 
R±i/L. 

It is convenient to proceed with the norm calculation 
in the selected class of functions not for the operator it- 
self R = Q^^^U but rather for its transposed counterpart, 
where operators Q'^^^ SloAU are rearranged. Considering 
that the norm squared of the chosen trial function is of 
the order of i, we actually need to estimate (by an order 
of magnitude) the expression 



|R|P = sup - dz y] 
qeM+i/L L Jc . 

X G 



^""^ (z,zi)G$.^'(z,Z2)e-"^(^^-^^) 



dzidz2U*j^{z)Ukniz) 
(A2) 



The potential AVn{z), see Eq. ( |l9b[ ), does not mod- 
ulo exceed the average potential Vn{z). Hence, we can 
simplify the integrand in Eq. ( |A2|) b y setting all the fac- 
tors g{z\n) = 1 in expression (|30|) and replacing the 
quantity >f„ in the exponentials by the wavenumber 



kn ~ Y " Vn{z). The reflection coefficient (32a) is not 
modulo larger than unity (see Fig. ^ , so instead of exact 
expression ( pO| ) we will use in what follows its truncated 
version, namely. 



■ exp 



(A3) 



With this expression, only two first terms in Eq. ( |30| ) 
are actually allowed for, whose amplitude factors are 
q["^ and ^22^- addition, in the right-hand sides of 



Eqs. (31a) and (31b) we retain the terms without the fac- 
tors of TZn, bearing in mind that the order-of- magnitude 
estimation only is of importance. All the unrecorded 
terms contribute to the norm (Al) additively with func- 
tion (A3), so they cannot affect significantly the esti- 
mates obtained below. 



As far as the intermodc potentials in Eq. ( |A2D are con- 
cerned, their expression ( |6b| ) contains two functionally 
different terms (with zero and non-zero average value). 
We will consider their contribution to the mode-mixing 
operator norm separately. The contribution of the first, 
quadratic in w' (z) term, after the function ( [A3| ) being 
represented in terms of the Fourier integral, is written as 



ml 



sup ^ ( — ) iQnl^ iBinDin){BknDkn) 



dsxdsi 



dz- 



-i{si—S2)z 



(*1 ^n*)(*2 ^n) 
I 



dzidz2e-*(«-'*i)^ie*(«-''2)^= . (A4) 



Here, the integral over z is calculated without difficulty 
if the length L is sufficiently large. Specifically, when 
\kn\L ^ 1 this integral may be thought of as the prelimit 
(5-function, 



A(A^) 



dz exp(±iK2:) ~ 



sin(KL/2) 



27r(5(K) 



Computing then the remaining integrals in Eq. (A4) 
and taking into account that the major contribution to 
the sums over i and k comes from the terms with the 
extended-type trial Green functions, we arrive at the es- 



timate 



IR 



■~4 



L 

Wo 



1 + S2 



(kLf . (A6) 



The contribution to norm (A2) of the second, linear 
in w'{z) term of potential ( ]6b[ ) is calculated in a similar 
manner, and the estimate of this contribution reads 



R| 



L 

Wo 



"2 



1 + s: 



■ikLf 



L 

^ (A7) 

By comparing the right-hand sides of Eqs. (A7) and ( [Aq ) 
one can reveal that their ratio equals to the value of 
(1 -I- S^)/(Sfc(i)^. Under the conditions suggested by in- 
equalities Eq. (E9), this parameter is greater than unity. 
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so the result (A7) may be reckoned as the final estimate 
of the operator R norm in the slope parameter region 
S < Scr, where the intermode scattering is ensured by 



extended modes. In the region S > Scr all the trial 
modes become evanescent, so it does not make any sense 
to take account of their mutual intermode scattering. 
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